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Based on an improved isospin- and momentum-dependent interaction, we have studied the quali¬ 
tative effect of isospin splitting of nucleon effective mass on the specific shear viscosity of neutron-rich 
nuclear matter from a relaxation time approach. It is seen that for m * > m p , the relaxation time of 
neutrons is smaller and the neutron flux between flow layers is weaker, leading to a smaller specific 
shear viscosity of neutron-rich matter compared to the case for m„ < m*. The effect is larger in nu¬ 
clear matter at higher densities, lower temperatures, and larger isospin asymmetries, but it doesn’t 
affect much the behavior of the specific shear viscosity near nuclear liquid-gas phase transition. 
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Understanding the basic strong interaction and the 
properties of nuclear matter is the main purpose of nu¬ 
clear physics. The knowledge of transport properties of 
the hot dense matter is important in understanding the 
dynamics in heavy-ion collision experiments as well as 
the properties of proto-neutron stars. The quark-gluon 
plasma (QGP) produced in ultra-relativistic heavy-ion 
collisions is believed to be a nearly ideal fluid and has a 
very small specific shear viscosity 77/s [lrBl, he., the ratio 
of the shear viscosity 77 to the entropy density s. It has 
been further found that the 77/s decreases with increas¬ 
ing temperature in the hadronic phase while increases 
with increasing temperature in QGP, resulting in a min¬ 
imum value at the temperature of hadron-quark phase 
transition @, 1 - At even lower temperatures, the 77/s of 
nuclear matter with nucleon degree of freedom has been 
investigated from the relaxation time approach 01 and 
transport model studies fToi [bl . Similar to the behavior 
near hadron-quark phase transition, the 77/s also shows a 
minimum in the vicinity of nuclear liquid-gas phase tran¬ 
sition from various approaches [III- 1 I 6 I]. Since the cor¬ 
relation between the elliptic flow and the specific shear 
viscosity seems to be a general feature in both relativis¬ 
tic 0 and intermediate-energy heavy-ion collisions 0 , 
in the future it might be promising to measure the 77/s ex¬ 
perimentally, meanwhile providing an alternative way of 
searching for nuclear liquid-gas phase transition in heavy- 
ion experiments. 

In our previous studies, the specific shear viscosity of 
neutron-rich matter was investigated based on an isospin- 
and momentum-dependent interaction BEL. Recently, 
this interaction has been further improved IT 8t by intro¬ 
ducing more parameters so that detailed isovector prop¬ 
erties can be studied more flexibly. One of the isovec¬ 
tor properties is the neutron-proton effective mass split¬ 
ting, and it becomes recently a hot topic. The interest 
was inspired by the recent experimental data of double 
neutron/proton ratio from the National Superconduct¬ 
ing Cyclotron Laboratory, which seems to favor a smaller 


neutron effective mass than proton based on the calcu¬ 
lation using an improved quantum molecular dynamics 
model [l9l ]. However, the well-known Lane potential, rep¬ 
resenting the nuclear symmetry potential, i.e., the differ¬ 
ence between the mean-field potential of neutrons and 
protons, decreases with increasing nucleon energy, lead¬ 
ing to a larger neutron effective mass than proton [201 ]. To 
explore the possible uncertainty of neutron-proton effec¬ 
tive mass splitting on the 77 /s of hot neutron-rich nuclear 
matter, we extend our study of the specific shear viscos¬ 
ity with the improved isospin- and momentum-dependent 
interaction (ImMDI) in this brief report. 

The single-nucleon mean-field potential of the ImMDI 
interaction is written as [l8j 
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where r = 1(—1) for neutrons (protons) is the isospin 
index, p n and p p are number densities of neutrons and 
protons, respectively, <5 = (p n — p p )/p is the isospin asym¬ 
metry, with p = p n + Pp being the total number density, 
and f T (p) is the phase-space distribution function. The 
x parameter is used to mimic the slope parameter of the 
symmetry energy at saturation density po , while addi¬ 
tional two parameters 77 and z are introduced to adjust 
the symmetry potential U sy m at infinitely large nucleon 
momentum and the value of symmetry energy E sym at 
saturation density, respectively, and they enter the func¬ 
tional through the relations 
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where pfo is the nucleon Fermi momentum in symmetric 
nuclear matter at saturation density. The values of pa¬ 
rameters A;o, A u o, i?, Cio, C u o, A, and a as well as the 
corresponding macroscopic quantities of nuclear matter 
from the ImMDI interaction can be found in Ref. [l8| . 



FIG. 1: (Color online) Binding energy in symmetric nu¬ 
clear matter (a), symmetry energy (b), symmetry potential 
at saturation density (c), and relative neutron-proton effec¬ 
tive mass splitting (d) in nuclear matter at saturation density 
and isospin asymmetry 5 = 0.5 from the two parameter sets 
based on the ImMDI interaction. 


The values of x, y, and z do not affect the isoscalar 
properties of nuclear matter, and the binding energy in 
symmetric nuclear matter is shown in panel (a) of Fig. [Q 
The parameter z is set to be 0 in the present study, while 
x and y change respectively the magnitude and the mo¬ 
mentum dependence of the symmetry potential, and they 
both contribute to the density dependence of the sym¬ 
metry energy. With different combinations of x and y , 
one can get very similar density dependence of symmetry 
energy but different momentum dependence of the sym¬ 
metry potential, or equivalently, the isospin splittings of 
nucleon effective mass, as can be seen from panels (b), 
(c), and (d) in Fig. [Q with the effective mass calculated 
from 
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The parameter sets [(a; = 0), (y = —115 MeV)] with 
m* > m* and \{x = 1), (y = 115 MeV)] with m* < m* 
are thus chosen in the following study. 

We now briefly review the main ingredient of the re¬ 
laxation time approach used in previous studies 1,0- 


The shear viscosity is calculated by assuming that in the 
uniform nuclear system there exists a static flow field in 
the 2 direction with flow gradient in the x direction. The 
shear force, which is related to the nucleon flux as well 
as the momentum exchange between flow layers, is pro¬ 
portional to the flow gradient, and the proportionality 
coefficient, i.e., the shear viscosity, turns out to be [§] 
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In the above, d = 2 is the spin degeneracy, p = 
\Jpx + Py + Pz is the nucleon momentum, and n T is the 
local momentum distribution 
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with p T and T being the chemical potential and the tem¬ 
perature, respectively. p x /m* is the nucleon velocity be¬ 
tween flow layers. t t (p) is the relaxation time for a nu¬ 
cleon with isospin r and momentum p, and it can be 
further expressed as 

1 - 1 1 
Trip) T° ame {p) + Tr iff {p)' 

with r * ame ( di //) ^ b e j n g average collision time for a 
nucleon with isospin r and momentum p when colliding 
with other nucleons of same (different) isospin. For more 
detailed derivations and the expressions of the relaxation 
time, we refer the readers to Ref. (j|. The relaxation 
time depends not only on the medium properties such as 
the density, temperature, and isospin asymmetry, but on 
the nucleon-nucleon scattering cross section as well. The 
free-space proton-proton and neutron-proton cross sec¬ 
tions (crjy^f) are taken as the parameterized forms from 
Ref. |2lf . while the in-medium cross section is modified 
by the effective mass through [22f 

( 10 ) 

where y NN ip* NN ) is the free-space (in-medium) reduced 
mass of colliding nucleons. The reduced mass scaling of 
the in-medium cross section comes from the fact that the 
differential cross section is inversely proportional to the 
relative velocity between the two colliding nucleons HI] , 
while the difference between the scattering T matrix in 
free space and in the nuclear medium is neglected in the 
present qualitative study. 

Figure [5] displays the total relaxation time and those 
for nucleon scatterings with same or different isospins. As 
expected, all relaxation times are larger with in-medium 
cross sections compared to the results with those in free 
space. In neutron-rich nuclear matter, the scatterings 
are more frequent for neutron-neutron than for proton- 
proton, and protons have more chance to collide with nu¬ 
cleons of a different isospin than neutrons. The total neu¬ 
tron relaxation time dominates the shear viscosity due to 
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FIG. 2: (Color online) Momentum dependence of the total 
relaxation time and that for a nucleon to collide with other 
ones of same or different isospin in nuclear matter at satura¬ 
tion density and isospin asymmetry of <5 = 0.5 using free-space 
or in-medium nucleon-nucleon scattering cross sections. 


the sharper neutron momentum distribution in neutron- 
rich nuclear matter as can be seen from Eq. 0, and it 
is larger for m* < to* than for to* > m* as a result of 
the isospin-dependent modification from the in-medium 
effective mass, while the difference of the total neutron 
relaxation time for different isospin effective mass split¬ 
tings is smaller with free-space cross sections. 



10 20 30 40 50 1.0 1.5 2.0 

T (MeV) p/p 0 


FIG. 3: (Color online) Specific shear viscosity with free-space 
((a), (b)) and in-medium ((c), (d)) cross sections in nuclear 
matter of isospin asymmetry 5 = 0.5 at different densities and 
temperatures for m* > m* and < rrip. 

The extensive results of the specific shear viscosity in 
nuclear matter of isospin asymmetry S = 0.5 at vari¬ 


ous densities and temperatures are shown in Fig. [3j Al¬ 
though the momentum occupation probability n T (p) at 
finite temperature depends on the neutron-proton effec¬ 
tive mass splitting [18|] , the difference of the entropy den¬ 
sity, which is calculated from 

f 0 (3 

s= J [n T In n T + (1 - n T ) ln(l - n T )\ , (11) 

turns out to be small between to* > to* and to* < to*. 
Despite of the similar relaxation time for different isospin 
splittings of nucleon effective mass using free-space cross 
sections, a larger neutron effective mass leads to smaller 
neutron fluxes between flow layers, reducing the shear 
viscosity as can be seen from Eq. 0. We note that this 
is a robust feature even if we use a naive mean-free-path 
formular. The isospin-dependent modification for the in¬ 
medium cross sections further enhances the difference of 
the specific shear viscosity between to* > m* and to* < 
to* by giving a smaller relaxation time for neutrons in 
the former case, as discussed in Fig. [2] In addition, one 
sees that the difference is larger at higher densities and 
lower temperatures when the relative isospin splitting of 
nucleon effective mass is generally stronger |la |. 




FIG. 4: (Color online) Isospin asymmetry dependence of the 
specific shear viscosity in nuclear matter at different temper¬ 
atures and densities for m > rrip and m), < rrip. 

The dependence of the specific shear viscosity on the 
isospin asymmetry is displayed in Fig. U at different tem¬ 
peratures and densities. It is seen that 77 /s increases with 
increasing isospin asymmetry 5 faster than a parabolic 
relation, and it increases even faster at higher densities 
or lower temperatures. One thus expects that the effect 
discussed here is quite relevant for the evolution of hot 
neutron stars with a large neutron excess. 

During the liquid-gas phase transition (LGPT) in nu¬ 
clear matter, each phase satisfying the Gibbs condi- 
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FIG. 5: (Color online) Temperature evolution of the entropy 
(upper panels) and specific shear viscosity (lower panels) in 
the presence of nuclear liquid-gas phase transition (LGPT) at 
fixed external pressure P = 0.05 ((a), (d)), 0.10 ((b), (e)), 
and 0.15 MeV/fm 3 ((c), (f)) and isospin asymmetry 5 = 0.5 
for rrin > m * and m * < rrip. 

tion [24j, [ 2 ^] has its own volume fraction. The thin lines 
in panels (a), (b), and (c) of Fig. [5] show the evolution 
of the entropy when the nuclear matter of isospin asym¬ 
metry S = 0.5 is heated at different fixed pressures. If 
the occurrence of the nuclear LGPT is taken into ac¬ 
count, the entropy evolution will follow the thick lines, 
with the overall entropy density from that in each phase 
weighted by the volume fraction. In infinite nuclear mat¬ 
ter the total shear viscosity can also be calculated from 
that in each phase weighted by the volume fraction [li| , 
and the temperature evolution of the specific shear vis¬ 
cosity is shown in panels (d), (e), and (f) of Fig. [5] A 
minimum value of rj/s is seen at a higher temperature 
with increasing pressure, and this value is also smaller at 
a larger external pressure. The isospin splitting of nu¬ 


cleon effective mass has very small effects on the entropy 
evolution, while the specific shear viscosity is smaller for 
m* > to* than for to* < to* at lower temperatures 
(higher densities) in the liquid phase side, but the dif¬ 
ference is negligible at higher temperatures (lower densi¬ 
ties) in the gas phase side. Moreover, the minimum point 
of r)/s is not affected by the isospin splitting of the nu¬ 
cleon effective mass. This general feature, which is not 
sensitive to detailed nuclear interaction, might be help¬ 
ful in searching for the occurrence of the nuclear LGPT 
in low- and intermediate-energy heavy-ion collisions if 
people find ways to measure the specific shear viscosity 
there [l7| . 

In summary, the specific shear viscosity with differ¬ 
ent isospin splittings of nucleon effective mass has been 
studied in neutron-rich nuclear matter based on an im¬ 
proved isospin- and momentum-dependent interaction. 
Qualitatively, it is seen that the specific shear viscos¬ 
ity is larger for to* < to* than for to* > m*, and 
the difference is more obvious at higher densities, lower 
temperatures, and larger isospin asymmetries. This is 
due to different neutron fluxes between flow layers as 
well as the isospin-dependent modification to the in¬ 
medium nucleon-nucleon cross sections. On the other 
hand, the behavior of the specific shear viscosity near 
nuclear liquid-gas phase transition remains robust and 
seems to be insensitive to the detailed nuclear interac¬ 
tion. Our study may be helpful in understanding the 
transport properties of the hot neutron-rich nuclear mat¬ 
ter produced in heavy-ion collision experiments as well as 
that in hot neutron stars. 
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